Exact mean first-passage time on the T-graph 
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We consider a simple random walk on the T-fractal and we calculate the exact mean time r 9 
to first reach the central node io. The mean is performed over the set of possible walks from a 
given origin and over the set of starting points uniformly distributed throughout the sites of the 
graph, except io. By means of analytic techniques based on decimation procedures, we find the 
explicit expression for r 9 as a function of the generation g and of the volume V of the underlying 
fractal. Our results agree with the asymptotic ones already known for diffusion on the T-fractal and, 
more generally, they are consistent with the standard laws describing diffusion on low-dimensional 
structures. 

PACS numbers: 05.40.Fb 



I. INTRODUCTION 

Many problems in physics and chemistry are related to 
random walks on fractal structures [jj, [2|, |3| • The main 
reason is that such structures are able to mimic the inho- 
mogeneity and scale-invariance typical of disordered ma- 
terials. An important class of fractal structures is given 
by the so called exactly decimable fractals which include 
deterministic finitely ramified fractals such as the Sier- 
pinski gasket and the T-graph. These structures, being 
amenable to renormalization procedures, allow exact an- 
alytic calculations @, H, [1] . 

In general, the lack of translational invariance implies 
significant corrections to the standard laws describing dif- 
fusion on regular lattices. Indeed, in the fields of reaction- 
diffusion and transport theory, a question of longstand- 
ing interest concerns the interplay between spatial extent 
and system dimensionality in affecting the reaction kinet- 
ics and the transport efficiency Q. 

A fundamental quantity characterizing diffusion is the 
mean first-passage time (MFPT), i.e. the expected time 
for a random walker, starting with equal probability at 
any site i ^ io, to first reach a given site io- This prob- 
lem was first set up by Montroll [8| in the case of regular 
structures and later extended to more complex substrates 
B 03, HO, H Gl EI HI- Notice that this definition of 
MFPT involves a double average: the first one is over all 
the walks from a given origin i, then you must average 
over a uniform distribution of initial sites, whose support 
is the whole set of graph sites, except io- 
The mean first-passage time is also intrinsically related 
to a number of different problems [HJ . In the context of 
reaction-diffusion processes it represents the mean time 
to react for a particle diffusing in the presence of an ac- 
tive site located in io [1, [H, [13, [H| , which is sometimes 
referred to as target problem. Not only: the mean first- 
passage time defined above also describes the asymptotic 
behaviour of the average time for two diffusive particles 
to first encounter [l9| . 

The MFPT has been previously studied on different 
kinds of structures and several analytical results have 
been found. Most of them consists of scaling relations 



FIG. 1: T-fractal of generation 4: V = 3 4 + 1. The next 
generation is obtained by performing the operation illustrated 
in the bottom on each bond. 



and asymptotic behaviours J^, [l3l [T3 |. while a very few 
exact results are known [Tl|,[l2]. Exact solutions on finite 
structures are especially longed for since they prove useful 
for a deeper comprehension of theoretical models and for 
checking approximate solutions or numerical simulations. 

Here, we derive the exact mean first-passage time for 
a random walker on the T-fractal, following the decima- 
tion procedure recently introduced by Kozak et al. and 
applied to the Sierpinski gasket [llj . In particular, we as- 
sume a simple random walker (RW) in the presence of a 
perfect trap fixed at the central site and we calculate the 
mean- walk length before absorption. The closed-form ex- 
pression we obtain for the latter, is akin to the one in II ill 
and consistent with known asymptotic results 0, ITol . Il9| |. 
Though they are both deterministic fractals, the Sier- 
pinski gasket and the T-fractal display significant dif- 
ferences: while the former models self-similar structures 
endowed with closed loops, the T-fractal is representative 
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for tree-like structures and bundled structures in general 
d, HE H3, HI] - For this reason it is worth extend- 

ing and comparing the related results. 

The paper is organized as follows. In Sec. |TT] we de- 
scribe the main features of the T-graph and we resume 
the analytic background underlying the analytic solu- 
tion; in Sec. IIIII we describe the decimation procedure 
applied and we obtain the exact formula for the mean 
first-passage time as a function of both the generation 
and the volume of the structure; finally Sec. IIVI includes 
conclusions and comments. 



II. THEORY 
A. Exactly Decimable Fractals 

A generic graph Q is mathematically specified by the 
pair {A, F} consisting of a non-empty, countable set of 
points, A joined pairwise by a set of links T. The cardi- 
nality of A is given by |A| = V representing the number 
of sites making up the graph, i.e. its volume. From 
an algebraic point of view, a graph Q = {A, T} is com- 
pletely described by its adjacency matrix A. Every entry 
of this off-diagonal, symmetric matrix corresponds to a 
pair of sites, and it equals one if and only if this cou- 
ple is joined by a link, otherwise it is zero. The number 
of nearest-neighbours of the generic site i, referred to as 
coordination number, can be recovered as a sum of adja- 
cency matrix elements: z% = J2jeA^v- These are used 
to build up the diagonal matrix Z%j = @. 
A very special class of graphs is given by the so called 
exactly decimable fractals which are geometrically invari- 
ant under site decimation. In general, all deterministic, 
finitely ramified fractals are exactly decimable. The solu- 
tion of both the random walk and harmonic oscillations 
problems can be obtained by standard renormalization 
group calculations based on real space decimation pro- 
cedures 0, H, @. The Sierpinski gasket, the T-graph 
(Fig. [T|) , the branched Koch curves are examples of ex- 
actly decimable fractals, which accounts for their popu- 
larity. Notice that such structures are characterized by 
strong restrictions on their topology which c an g ive rise 
to properties far from holding for all fractals [2(| ■ 

Here, we consider a T-fractal which is iteratively con- 
structed by performing the operation illustrated in Fig. [T] 
on each link. The number of iterations is called the gen- 
eration g of the fractal. At the 5-th generation the cardi- 
nality of the set of nodes V(g), hereafter called volume, 
is given by: V(g) = \A g \ = 3 9 + 1. The T-fractal has 
fractal dimension df — ss 1.584 and spectral dimen- 
sion d — w 1.226. We recall that the former gives 
the dependence of the volume of the system on its linear 
size L: 

V(g)~lL(g)] d ' = (2°) d f, 
while the latter governs (among other phenomena) the 
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FIG. 2: (Color on line) T-graph of generation 2 (top) and 3 
(bottom) with volume V — 3 2 + 1 and V = 3 3 + 1, respec- 
tively. Each site has been labelled according to the procedure 
described in the text. Dotted links belong to F 9 and colored 
sites to A 9 ; more precisely, pink sites are in A g ext and blue 



long-time properties of diffusion on the graph. Indeed, if 
we consider a random walker starting from a given site 
i of the graph, the probability Pa(t) of returning back 
to the starting point at time t, at long times, follows the 
law 

P lt (t)~t- d / 2 . 

When d < 2 the random walker is said to perform a 
"compact exploration" of the space [24j since the fractal 
dimension of its trajectory is greater than the dimension 
df of the substrate. 

It is worth underlining that the T-fractal is irregular, 
that is the coordination number is site dependent. We 
can distinguish among "internal site" with coordination 
number z = 3 and "external site" with z = 1. We call 
Ai n t and A ex t the set of internal and external sites re- 
spectively. Obviously A 9 = A 9 int U A 9 ext . 

Figure [2] shows generations g — 2,3 and the labelling 
method adopted, which will be useful in the next sec- 
tion. In each new generation, we label only the new sites 
while old sites keep their own labels. Hence, at gen- 
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eration g + 1 we name new sites progressively, starting 
from V(g) + 1. The first sites to be labelled are the 
innermost, the last ones are those farthest from the cen- 
tral site io = 1. At each generation we can distinguish 
sets of equivalent sites which are labelled anti-clock wise. 
Due to the symmetry of the T-graph, the cardinality of 
such sets is always a multiple of 3; for example we have 
{2,3,4,8,9,10} g=2 ,{14,15,16} ff=3 . 

Before proceeding further, let us resume some facts 
concerning unbiased diffusion on a T-fractal of generation 
g represented by the adjacency matrix A g (henceforth we 
will omit the subscript g). We consider a simple RW, 
starting, at t = 0, from site i\ at each time step (taken 
to be unity) the particle jumps with equal probability 
to any of its nearest-neighbour sites. Being Pji{t) the 
probability of going from i to j in t steps, the following 
Master equation holds: 

V 

^(i+l) = E(^ 1 ^ P ^W' (!) 

k=l 

which states that at each time step, the jumping proba- 
bility from an internal site is ^ . From the previous equa- 
tion it follows that Pji(t) — [(Z" 1 A) t ]ji. It is also easy 
to verify that the Markov chain representing such a ran- 
dom walk is ergodic and the particle will visit all sites 
with probability 1, independently of the origin i. Con- 
sequently, the walker will reach any site with probability 
1, in a time possibly diverging when g — > oo. 

B. Mean Time to Absorption 

Let us consider a perfectly absorbing trap, fixed on the 
central site (labelled with index io = 1) of the T-graph. 
Our aim is to obtain, through a decimation procedure, 
an exact, closed- form solution for the average time to 
absorption, where the average is meant both over all the 
possible walks starting from the same origin i and over 
all sites i ^ io taken as origin of the walk. 
The special choice we made for the trap location makes 
the decimation procedure easier to be applied as we can 
identify the site io since the first generation. 

We now introduce t? q to be the q-th (q — 0, 1, ...) mo- 
ment of the trapping time for a walk starting from i on 
the g-th. generation of the graph. Obviously, regardless 
of g, Tq — and t 9 = 1. The latter is the zero-th mo- 
ment of the distribution of the time to absorption given 
i as origin and it is unitary because the walker will be 
trapped with probability 1, whatever its origin. 

The starting point for our analytic treatment is the 
discrete differential equation introduced in pj| : 

V(g) 

-E A ^+i = (9 + iK,, (2) 

3=2 

where A = AZ~ Y — I is a normalized version of the dis- 
crete Laplacian whose first row and column (correspond- 



ing to the trap site) have been removed. We recall that A 
is a nonsingular matrix and each row has sum zero, apart 
from those corresponding to sites V(g—1) + 1, V(g — 1) + 2 
and V(g — 1) + 3, i.e. the three nearest-neighbours of the 
trap, for which the sum is — ^. As discussed in [7|, Eq. [2] 
can be generalized to the case of two or more particles 
simultaneously diffusing. 

In the following, we just focus on the set of first mo- 
ments Tji, for which Eq. |(5J) simplifies into 

EA- T 9 - T 9 - 1 

3 

Henceforth we can drop the index corresponding to q 
without ambiguity: 

-£A y rf = l. 

3 

Now, we implement the average over the starting site i, 
chosen according to a uniform distribution in A 9 \ {io}' 

. Via) V{g)V(g) 

VJ > i=2 yi " i=2 j=2 

(3) 

In the next Section we derive some recurrence relations 
which allow to simplify the previous equation. 

Equation [3] can be very easily interpreted if we look 
at the random walk as a Markov chain. Indeed, — A -1 
is just the fundamental matrix for the process, whose 
entry i,j represents, by definition, the expected number 
of times that the process is in the transient state j, being 
started in the transient state i. 

Finally, notice that, due to the symmetry and the ab- 
sence of loops characterizing the structure under con- 
sideration, the mean time to absorption found in this 
case just corresponds to the mean time to reach either 
site 2, 3, 4 on a T-fractal of generation g — 1. Otherwise 
stated, if we call it 9 the mean time to first reach site 2 
(or, analogously, 3,4), then it 9 = t 9+1 . 

III. DECIMATION PROCEDURE 

The number of terms to sum up in Eq. [3] grows expo- 
nentially with g, hence, a direct calculation of rf and t 9 
can be accomplished straightforwardly only for the very 
first generations (see Tab. HI and HI)) . Such data allow to 
get some recurrence relations useful for the derivation of 
the final formula. First of all, notice that for a given site 
i, we have t? + = 6t 9 : in each generation the chemical 
distance from i to the trap doubles while the mean time 
to first reach the trap increases by a factor 6. This ex- 
act scaling follows from the symmetry and decimability 
of the graph and it is consistent with the random walk 
dimension on the T-fractal: d w = = ^5 fij]. 

_ d _ log 2 ! 

Furthermore, at each generation g we insert on the ex- 
isting fractal some new links and some new vertices. We 
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TABLE I: Mean time to absorption t? for a random walker starting from a given site i. For these values the average is only 
performed over all possible random walks sharing the same origin. Due to the symmetry of the T-fractal we can distinguish 
sets of equivalent sites such that, if taken as origin of the walk, they provide the same mean time rf . Notice that, for the 
farthest sites from the trap, having chemical distance 2 9_1 , the average absorption time is 6 9_1 , consistently with the random 
walk dimension d w . 

g\j (2, 3, 4) (5, 6, 7) (8, 9, 10) (11, 12, 13) (14, 15, 16) (17, 18, 19, 20, 21, 22) (23, 24, 25, 26, 27, 28) (29,30,31) (32, 33, 34) 
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210 


216 


53 


54 



call such sets A 9 = A^A 9 " 1 and V = r 9 \r 9 -\ respec- 
tively. It is easy to see that |A 9 | = V(g) — V{g — 1) = 
2 • 3 9 " 1 and 2 • |f 9 | = |A 9 |. Moreover, for each new link 
added we have a new couple of connected vertices j e xt 
and ji n t , belonging to A 9 ext and A 9 nt , whose coordination 



numbers are 



1 and Zj int = 3, respectively and: 



A 9 



A 9 1 1 \ 9 



-A 5 = A g ext = Af nt . 



(4) 



For example, A e 



{14, 15, 16, 23, 24, 25, 26, 27, 28}, as 



shown in Fig. [51 Now, it is easy to see that 

1 (5) 



T? = 7-? 



since a RW starting from i ext is necessarily on i int at 
time t = 1. These facts hold regardless of the generation 
g and for any connected couple chosen from A 9 . Thus we 
can write 



V(g) 

E 

i=2 



E 

ieAs- 1 



(6) 



ieAs 



In the last sum we can highlight the contribution from 
external and internal sites and, exploiting Eqs. 2] and El 



E-f 



E 



^^-2^ rf + |AL|. (7) 



iGA 9 



Let us now focus on the sum appearing in the right-most 
side and estimate it in the case g = 3 depicted in Fig. [5] 
The mean time to absorption for a RW starting from site 
5 can be expressed as: 



(1 



(1 



(1 



(8) 



'14 "r '23 



'24 



where Pfci represents the transition probability from state 
j to state k. Hence, rf is just the average of the absorp- 
tion times from 11, 17, 18, which mirrors the barycentric 



position of 5 with respect to the latters. Also notice 
that site 5 belongs to Af~ t . Since these facts hold for 
any analogous subtree of generation 2, we are allowed to 
write: 



E ^ = 3 - 6 E 



From Eqs. [7] and [9j 
E rf 



ieA" 



-3 9 - 1 = is j2 r r 1 . 



(9) 



(10) 



and solving this recurrence relation we obtain: 

33-1 



E 



-(1 + 4-6 9 " 1 ). 



(11) 



By plugging the last expression into Eq. ([6]) we get 
(V(g) - l)r 9 =6(V(g - 1) - l)^ 1 -!- 

+ _ .39-1(1 +4. 6 9 - 1 ) + 3 9 - 1 
5 



and, dividing by V(g) — 1 = 3 9 : 

1 



2r 9 ~ 



15 



(7 + 8 -6 9 " 1 ). 



(12) 



The last expression is, again, a recursive equation, whose 
solution provides the exact time to absorption: 



4(-7 + 5-2 9 + 2-6 9 ). 
15 



(13) 



A numerical check of this formula can be attained by 
comparing r 9 , 1 < g < 6 with data obtained by direct 
calculation and reported in Tab. |TTJ the agreement is 
perfect. 

It is also possible to obtain an expression for r 9 as a 
function of the volume V(g) = 3 9 + 1. In fact, recalling 
the spectral dimension for the T-graph, we can write 2 9 = 
(V(g) - 1) 2 /<H and 
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5{V(g) - l) 2 ^" 1 + 2{V{g) - l) 2/i - 7 



(14) 
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TABLE II: Mean first-time r 9 obtained by direct calculation 
from Eq. [3] since this implies a sum over a number of terms 
exponentially increasing with g, only small generations have 
been considered. 

1 Vfe) ~? 



4 
10 
28 
82 
244 
730 



51/9 
837/27 
14391/81 
254421/243 
4550175/729 



Notice that the last expression gives the exact, explicit 
dependence of t 9 on V(g). In the asymptotic limit, as 
V(g) diverges, 



V(g) 



2/d 



(15) 



This result is consistent with the leading behaviour of t 9 
on the T-fractal discussed in [9j and, more generally, to 
the leading behaviour of the trapping time on low dimen- 
sional (d < 2) structures. In fact, as already remarked, 
the MFPT represents the mean trapping time rt rap for a 
diffusing particle in the presence of a fixed perfect trap 
(or, symmetrically, the trapping time for an immobile 
target in the presence of a diffusive trap). Being P s (t) 
the survival probability at the £-th step P s (t), i.e. the 
probability that the RW has not yet reached the trap 
site, then [![: 



Ttrap 



dP s (t) 
dt 



t dt. 



The asymptotic expression for the survival probability is 



with 



t d ' 2 , d<2 
t, d>2 



(16) 



(17) 



One therefore expects that, for low dimensional struc- 
tures, Ttrap ~ V 2 / d . Our result and the one in [ll[ estab- 
lish this relationship rigorously. 

As can be shown by scaling arguments, on long times, 
the factor 6 9 ~ V(g) 2 ^ d , which is the leading term in 
Eqs.[l3l[T4l is involved in all the dynamical properties of 
diffusion on the T-fractal. For example, the characteristic 
time in the exponential decay of the survival probability 
increases with the generation of the tree as 6 9 @ . 

Equation [14] is also consistent with a recent result ob- 
tained by Condamin et al. [l5| who found the asymptotic 
(large V) expression for the average time r r taken by a 
RW on a generic scale-invariant structure to first reach a 
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FIG. 3: (Color on line) Mean first-passage time for a simple 
random walker moving on a T-fractal (black stars) and on 
a Sierpinski gasket (red circles) as a function of the volume 
Vspan (the volume the RW can actually span before being 
trapped), according to exact analytic solutions. 



trap distant r from the starting point. In particular, for 
compact exploration (d w > df), r r ~ V ■ r du >~ df . Thus, 
if we fix the trap on a particular site io, we can obtain 
an estimate for the MFPT, by simply averaging over all 
possible distances r from io- For the T-fractal considered 
here, io = 1 and we can write: 



V(g) 

where n(r) is the number of sites distant r from io and 
L(g — 1) is the largest distance from the central site. 
Under the above mentioned assumption of large volume, 
we can adopt a continuous picture and n(r) ~ 3-r rf/_1 ; by 
integrating the previous expression, we get t 9 ~ V(g) 2 ^ d , 
as expected. 

Finally, in Fig. [3] we compare Eq.[TT]with the analo gou s 
formula found for the Sierpinski gasket (d = p^-t) in pl| : 



r 9 - 



(2V{g)-3)l , 2(2V(g)-3)< 



In the asymptotic limit, t 9 diverges faster for the T- 
fractal. This can be trivially drawn algebraically while, 
from a topological point of view, it evidences the role of 
loops in reducing the average distance between two ran- 
dom sites, making the diffusive particles survive shorter. 
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IV. CONCLUSIONS 

In this work we study the mean first-passage time t 9 
for a random walker on a T-fractal. The latter, being ex- 
actly decimable, allows the use of the powerful technique 
of exact renormalization. We find an exact, closed-form 
solution for t 9 as a function of either the generation g 
and the volume V(g). The leading term of t 9 is consis- 
tent with known asymptotic results. It should be under- 
lined that an exact solution on a finite system is generally 
very useful in order to understand more quantitatively 
the asymptotic limit. 

Our findings are interesting also in the light of a recent 
result concerning the survival probability for the trapping 



problem A + B — ► B, with both species, A and B, dif- 
fusing There, it was shown that on low-dimensional 
structures (d < 2) the survival probability for an A par- 
ticle asymptotically does not depend on its diffusivity 
constant Da- Otherwise stated, at long times, the tar- 
get problem and the trapping problem provides the same 
results. Hence, exact results concerning the target prob- 
lem, also provide the correct asymptotic behaviour for 
the trapping problem with diffusive traps. 



The author is grateful to D. Cassi and R. Burioni for 
useful discussions and comments. 
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